Entropic attraction of adhesion bonds toward cell boundaries 



O 

(N 

> 
O 

ON 



o 

C3 



I 

C 

o 
o 



(N 
> 
00 

in 
in 

d 



X 



Noam Weil and Oded Farago 1 
'Department of Biomedical Engineering, Ben Gurion University, Be'er Sheva 84105, Israel 

Adhesion bonds between membranes and surfaces are attracted to each other via effective in- 
teractions whose origin the entropy loss due to the reduction in the amplitude of the membrane 
thermal fluctuations in the vicinity of the adhesion bonds. These fluctuation-induced interactions 
are also expected to drive the adhesion bonds toward the rim of the cell, as well as toward the 
surfaces of membrane inclusions. In this paper, we analyze the attraction of adhesion bonds to 
the cell inner and outer boundaries. Our analysis shows that the probability distribution function 
of a single (diffusing) adhesion bond decay algebraically with the distance from the boundaries. 
Upon increasing the concentration of the adhesion bonds, the attraction to the boundaries becomes 
strongly self- screened. 

PACS numbers: 



I. INTRODUCTION 

Biological membranes are essential to every living or- 
ganism. They consist of a bilayer of amphiphilic lipids, 
proteins and cholesterol, and serve as a barrier between 
the cytosol and the extracellular matrix in biological cells 
In eukaryotic cells, membranes surround the or- 
ganelles, allowing for organization of biological processes 
through compartmentalization. Membranes are involved 
in a variety of cellular processes, including the transition 
of nutrients in and out of the cell, signal transduction 
and communication [2]. 

One of the processes taking place at the cell mem- 
brane is cell adhesion. This process, during which the 
cell membrane is attracted to another interface (which 
may be the membrane of another cell) can, in principle, 
be facilitated by non-specific attractive interactions (e.g., 
Coulomb and van der Waals interactions). Most often, 
however, cell adhesion is mediated by specific binding be- 
tween receptors that reside on the membrane surface and 
ligands located on the opposite surface [1] . Specific adhe- 
sion usually occurs in regions with high density of recep- 
tors and ligands. When facing a surface with enough lig- 
ands, the receptors may cluster into highly concentrated 
adhesion domains to establish strong binding 4, 5] . Spe- 
cific bioadhesion occurs in a variety of cellular processes, 
including binding of white blood cells to pathogens @, 
binding and fusion of drug carrier liposome to target cells 
Q, cadherin- mediated adhesion of neighboring cells 0, 
focal adhesion of cells to the extracellular matrix Q , and 
cell signaling. 

Adhesion induced domain formation results from both 
the activity of the actin cytoskeleton, and from attractive 
interactions which are mediated by the membrane curva- 
ture elasticity • The membrane- mediated interactions 
give rise to two distinct, but not entirely unrelated, mech- 
anisms facilitating cluster formation. The first mecha- 
nism is the Casimir-like attraction between two bonds, 
which originates from the reduction in the amplitude of 
the membrane thermal fluctuations in the vicinity of the 
bonds and the associated entropy loss. The resulting in- 
crease in the free energy of the membrane is minimized 



when the two bonds come within close proximity of each 
other, in which case the area where the thermal fluctu- 
ations are suppressed is reduced [ll|. The second mech- 
anism is related to the tendency of new bonds to form 
next to existing ones. This "binding cooperativity" ef- 
fect stems from the fact that the ligand-receptor binding 
probability is highest when the separation between the 
ligand and the receptor on the opposite surfaces is close 
to the bond length - conditions which are likely to occur 
near already existing bonds 12]. In addition, since the 
formation of new bonds leads to an overall reduction in 
the membrane roughness, it also increases the binding 
probability and assists the formation of even more bonds 
1 1 31 ] - Generally speaking, the first mechanism is likely to 
contribute to the aggregation of irreversible bonds with 
large binding energies, while the second one is probably 
more relevant to the formation of metastable domains of 
weaker reversible bonds. 

In a previous publication we proposed a framework 
for studying the Casimir-like interactions between adhe- 
sion bonds [14]. Our analysis is based on the following 
idea: The Casimir-like interactions are induced by the 
membrane thermal fluctuations whose amplitude is sup- 
pressed due to the presence of adhesion bonds. To calcu- 
late the free energy cost of a given distribution of adhe- 
sion bonds, a free energy density is assigned to each unit 
area of the fluctuating membrane, and the free energy 
density is integrated over the entire membrane area. The 
free energy density reflects the extend by which the am- 
plitude of the fluctuations is locally reduced compared to 
the fluctuations of a free (non-adhering) membrane. The 
key point of the method is the assumption amplitude of 
thermal fluctuations depends on the distance ci m i n to the 
nearest adhesion bond. Accordingly, the free energy den- 
sity associated with each unit area of the membrane is 
given by (see derivation in [l4j]): 



V(r) = 



(1) 



2 



The total attachment free energy is given by 

(2) 

where the integration is carried over the entire membrane 
area accept for the immediate vicinities of the adhesion 
bonds. The sum on the r.h.s. of eq. [2]is the discrete ana- 
log of the integral, which is useful for lattice simulations 
and numerical integrations. In the sum, the variable Sj 
takes the value s, = 1 for the sites occupied by adhesion 
bonds and Si = for an empty sites that represent fluc- 
tuating unit areas of size I 2 . Only the latter contribute 
to the sum. 

We studied the model defined by eq. [5] by using both 
a mean field analysis and Monte Carlo simulations flij . 
Our analysis revealed (in agreement with other theoret- 
ical studies [TH, [l6j]) that the fluctuation-induced inter- 
actions are not strong enough to allow the formation of 
adhesion domains. However, usually the adhesion bonds 
also interact with each other via short range attractive 
hydrophobic and depletion forces [l7|; and when these 
residual interactions are introduced into the model, con- 
densation of adhesion domains becomes thermodynam- 
ically possible [3, EH- Our analysis showed that the 
fluctuation-induced Casimir-like interactions greatly re- 
duce the strength of the residual attractive interactions 
required for cluster formation by almost half an order of 
magnitude, to below the thermal energy ksT. 

In this work we take a step forward and study the for- 
mation of adhesion domains in cellular membranes which 
also include trans-membrane proteins. Since these inclu- 
sions significantly influence the membrane spectrum of 
thermal fluctuations [ijj], their presence is expected to 
modify the fluctuation-induced interaction between ad- 
hesion bonds and their aggregation behavior. In this con- 
text, it is worthwhile to mention that there is a great bulk 
of experimental [l9l - [2~fl ] and theoretical [22| - [25j works on 
the fluctuation-mediated potential between the inclusions 
themselves. Here, we tackle a different problem and ex- 
plore how a single inclusion affects the aggregation of 
adhesion bonds. Since the formation of adhesion cluster 
is driven by the tendency to localize the reduction of the 
thermal fluctuations in a restricted area, it is natural to 
expect that adhesion bonds would exhibit affinity to the 
surface of the inclusion. Physically, the effect of an inclu- 
sion on the thermal fluctuations should be similar to that 
of a cluster of adhesion bonds of the same size. However, 
as discussed above, without direct residual interactions 
the system is in the gas-phase and a cluster of adhesion 
bonds is not stable. Therefore, the presence of the inclu- 
sion is not equivalent to the formation of a cluster, which 
is only transient. 

The paper is organized as follows: In sec. [TTJ we exam- 
ine the attraction between a circular inclusion and a sin- 
gle adhesion bond. We derive the probability distribution 
of the adhesion bond from the center of the inclusion. We 
show that in contrast to the weak attraction between two 




FIG. 1: Schematic of the system under investigation, consist- 
ing of a circular membrane of radius L with a circular inclu- 
sion of radius ro at the center. The membrane is supported 
by a flat surface to which it is attached by a single adhesion 
bond located at f. 

adhesion bonds (see ref. [26]), the fluctuation- induced at- 
traction between the adhesion bond and a large inclusion 
is sufficiently strong to keep the adhesion bond in the 
vicinity of the inclusion. In sec. IIII1 we study the distri- 
bution and aggregation of many adhesion bonds in the 
presence of the inclusion. Using Monte Carlo simulations 
with different boundary conditions and different adhesion 
bonds densities, we demonstrate that the strong attrac- 
tive potential calculated in sec. |H] is screened and felt 
only by the adhesion bonds located near the surface of 
the inclusion. In sec. IIV1 we summarize and discuss our 
results. 



II. THE BOND-INCLUSION PAIR 
INTERACTION 

We consider the model system shown schematically in 
fig. [TJ consisting of a circular membrane of radius L, a 
circular inclusion of radius ro <C L fixed at the center of 
the membrane, and a single adhesion bond located at r. 
The total attachment free energy of the system is given 
by eq.® 

F(?,r ,L) = f L 1 -dP, (3) 

Jr +i 7rrf^ in (r',r,r ,L) 

where the integration is carried over the whole area of 
the membrane, excluding regions of microscopic size I 
near the inner (|r'| = ro) and outer (|r'| = L) bound- 
aries, and in the vicinity of the adhesion bond (r' = r). 
In eq. 03 <i m i n denotes the distance of a given point r' 
on the membrane to the nearest "obstacle" , which may 
be either the adhesion bond or the surface of the in- 
clusion. The assumption underlying eq. [3] is that the 
membrane separation at the inclusion is the same as the 
length of the adhesion bond. This feature is also likely to 
contribute to the binding cooperativity mechanism men- 
tioned above in section^ which we do not consider in this 
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FIG. 2: Data collapse of the attachment free energy for differ- 
ent values of \r\, ro, and L, as a function of the dimensionless 
distance x — (|F| — ro)/(L — ro). The data in (A) and (B) cor- 
responds to systems with open and closed outer boundaries, 
respectively. 



work. Had there been a length mismatch between the in- 
clusion and the adhesion bonds, it would have involved 
a bending elasticity energy cost that would weaken the 
fluctuation-induced attraction. In a binary mixtures of 
adhesion bonds, the inclusion is likely to have a higher 
affinity to the adhesion bonds with the smaller length 
mismatch, which may lead to phase segregation in the 
mixture [HlH- 

In what follows, we consider membranes with both 
"open" and "closed" outer boundaries. At a closed 
boundary, the membrane is attached to a frame, whereas 
at an open boundary, the membrane fluctuates freely. 
Thus, the former is also considered as an obstacle, while 
the latter is not. The open and closed boundary condi- 
tions represent two limiting cases. Cellular membranes 
are attached to the actin cytoskeleton located mainly at 
the periphery of the cell., which corresponds to an inter- 
mediate case between open and closed boundary condi- 
tions. 

The integral in eq. [3] has been evaluated numerically 
for different values of ro and L, and for ro < \r\ < L. 
Our calculations reveal that the attachment free energy 
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FIG. 3: The width A* of the shell around the inclusion, as a 
function of the radius of the inclusion ro. The dashed hori- 
zontal line corresponds to the solution of eq. [9j A* = 21, for 
asymptotically large ro. 



is very well approximated by the following expression 

(H-r ) 



F(r,r ,L) ~ 2k B T\n- 



I 



+ nfc 5 Tln (L - |r|) +C(r ,L), (4) 

with n = 1 for an open outer boundary, and n — 2 for a 
closed boundary. Introducing the dimensionless distance 
< x = (\f\— ro)/(L— ro) < 1, cq.|4]can be also expressed 
as 



F(x) = F(r,r ,L) - In 



L - Tp 
I 



2+n 



-C{r ,L) 



2k B Tln(x) + nk B T\n(l - x) 



(5) 



Our numerical results along with eq. [5] are plotted in 
figs. EJA) and EJJB) for an open and closed boundary, 
respectively. When the adhesion bond is located far 
away from the outer boundary, i.e. for x <C 1, the 
attachment free energy is dominated by the first term 
F ~ 2k B Thi(x), which has been previously derived as 
the pair potential of mean force between two adhesion 
bonds [26]. In the vicinity of an outer closed bound- 
ary ( x ~ 1): one gets F ~ 2fcsTln(l — x), which im- 
plies that the adhesion bond is attracted to the nearest 
point on the outer closed boundary as if this point is 
another adhesion bond. For an open boundary, the pref- 
actor of the second term in eq. [5] is reduced by half from 
n = 2 to n = 1. This can be rationalized by noticing 
that an open boundary is identical to the midplane be- 
tween two adhesion bonds in a twice larger system. The 
pair potential between the adhesion bond and its image 
is 2fcBTln[2(l-x)] = 2fc B Tln[(l-x)] + C, and this en- 
ergy is equally divided between the real and imaginary 
halves of the system. 

How strong is the attraction between the adhe- 
sion bond and the inner inclusion? From eq. 0] we 



4 



find that the probability distribution function P(f) ~ 
cxp [—F(r, 7*0, L)/ksT] of the adhesion bond is given by 



P{f) = 



1 



1 



Z{\r\-r o y{L-\r\y 
where the normalization factor is 



Z 



L-l 



2irrdr 



,+, (r-r ) 2 (L-r)« 



(6) 



(7) 



In the case of an open outer boundary (n = 1), the 
mean distance between the adhesion bond and the sur- 
face of the inner inclusion is given by 



(r) - r = A 



2L\n(^^)+r (r -L) 



21+rg-L 
l(L-r a -l) 



(8) 



which for L 3> 7- = |r| 3> ro, simplifies to A ~ L/2. The 
fact that the mean separation A grows linearly with the 
size of the membrane may be considered as an indication 
that the adhesion bond is not bound to the inclusion. 
However, a closer inspection of eq. [6] reveals that the ad- 
hesion bond spends most of its time near the boundary 
of the inner inclusion. To quantify this phenomenon, we 
introduce the length scale A*, which is the width of the 
shell around the inner inclusion where the probability to 
find the adhesion bond is 0.5. The length A* is deter- 
mined by solving the equation 



Pace (A*) 



i-o+A* 



27r|7 ; '|P(|7 ; 1)dr 



2Lln(Al)+r (ro-L)(^) 



r +l 



2L\n(^f^)+r (r -L) 



21+rg-L 
l(L-r -l) 



0.5. 



(9) 



Our results for A* are plotted in fig.[3]as a function of r 
and for L ^> vq. The results shows that A* is a monoton- 
ically decreasing function of ro. For very small inclusions 
of size 7*o ~ 2, the length A* 3> 2. This result is consis- 
tent with the observation of ref. (2(| that the fluctuation- 
mediated pair potential is too weak to bind two adhesion 
bonds to each other. For inclusions of size ro > 20/, the 



length A* < 52; and for asymptotically large inclusions 
(ro —>■ oo), A* -*> 22. In other words, although A ~ L/2, 
the adhesion bond is likely to be found within a thin shell 
around the surface of a sufficiently large inclusion. 

For a closed outer boundary, the radial distribution 
function of the adhesion bond 



g(r) = 27rrP(|r|) 



(L-r ) 3 



(10) 



(The expression appearing in braces in eq. [10] is a nor- 
malization constant that depends on ro and L). The ra- 
dial distribution function g(r), which is plotted in fig. [4] 
for r = 41 and L = 100L, increases rapidly as one 
approaches both the inner and outer boundaries of the 
membrane. However, unlike the probability density per 
unit area P(r) which is symmetric with respect to the 
mid-radius r rn = (L + t*o)/2, the radial distribution func- 
tion g(r) ~ rP(|r|) is much higher near the outer bound- 
ary than the inner one. This, of course, is directly re- 
lated to the fact that there is simply more membrane 
area near the outer boundary. In the example plotted in 
fig-EK^o — 42, L — 1002), the probabilities to find the ad- 



hesion bond within a shell of size 42 around the outer and 
inner boundaries are 72.5% and 4.5%, respectively. For 
a bigger membrane of size L = 10002 (which, assuming 
that 2 is of the order of a few nanometers, is a reason- 
able estimate for the size of the cell plasma membrane), 
these probabilities change to 79% and 0.5%, respectively. 
These numbers suggest that the adhesion bond is likely 
to be "adsorbed" at the outer boundary. This, however, 
is only part of the story. If the adhesion bond is located 
near the inner boundary, it is not going to easily escape 
to the outer rim. The inset on fig. [4] shows the poten- 
tial of mean force F g {r) = — fcsTln \g (r)], exhibiting a 
free energy barrier of ~ 4fc bT that blocks the migra- 
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FIG. 4: The radial distribution function g as a function of r, 
for ro = 47 and i = 100/. The inset shows the corresponding 
free energy, F g — —kBT\n[g(r)l]. 



tion of the adhesion bond from the inner boundary to 
the outer one. The barrier on the opposite direction is 
of ~ IksT . For L = 1000/, the free energy barriers in- 
crease to ~ 7ksT and ~ 12fceT for the escape from the 
inner and outer boundaries, respectively. Thus, as in the 
case with an open boundary discussed above, an adhe- 
sion bond located near the inner boundary may remain 
trapped there for relatively long times. 



III. DISTRIBUTION OF ADHESION BONDS 
AROUND AN INCLUSION 

The previous section dealt with the localization of a 
single adhesion bond near the boundaries of the system. 
Cellular adhesion, however, involves an ensemble of ad- 
hesion bonds which may cluster into adhesion domains. 
In a previous work we demonstrated that the for- 
mation of adhesion domains cannot be induced by the 
Casimir-like interactions alone, and that it requires the 
existence of residual attractive interactions between the 
adhesion bonds. Here, we wish to understand whether 
this picture is changed by the presence of a membrane 
inclusion that serve as a potential nucleation seed for an 
adhesion domain. One can envision a scenario where a 
thin shell of adhesion bonds is formed around the in- 
clusion, effectively increasing its size and promoting the 
recruitment of more adhesion bonds. The opposite sce- 
nario, in which the addition of adhesion bonds self-screen 
the attraction toward the inclusion, is also plausible. In 
order to investigate this issue, we performed Monte Carlo 
lattice simulations of the model defined by cq. [2] on cir- 
cular systems with ro = 4/ and L = 100/ (where / is the 
lattice spacing). Our simulations were conducted at rel- 
atively low concentrations of adhesion bonds, cj) — 0.01 
and <fi = 0.05, with both open and closed boundary con- 
ditions. 

We measured the number density (per unit area) of 
adhesion bonds as a function of r. To allow compar- 
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FIG. 5: The normalized number density p(r) as a function 
of r = \r\. Data shown in the figure corresponds to ro = 4/ 
and L = 100/, for systems with (A) open and (B) closed outer 
boundaries. The solid lines depict the probability distribution 
of a single adhesion bond. 



ison with the probability distribution P{r) (eq. [B]), we 
normalized the number density to the total number of 
adhesion bonds. Our results for the normalized num- 
ber density, p(r), are depicted in fig. [S]A) andJSfB) for 
open and closed boundaries, respectively. Both figures 
demonstrate that the fluctuation-induced attraction to 
the boundaries is self-screened by the presence of adhe- 
sion bonds. Comparison of our results for cj) = 0.01 (de- 
picted by stars in both figures) and <j> = 0.05 (pluses), 
with eq. [6] for a single adhesion bond (solid lines) , re- 
veals that the self-screening effect is increased with the 
increase in the concentration <j). For an open boundary, 
the slight increase in P(r) near the outer boundary disap- 
pears already at = 0.01. For <j> = 0.05, the density p{r) 
is almost constant with no significant affinity of adhesion 
bonds toward the inner boundary. The same trends are 
also observed in the case of a closed boundary, where 
P{r) increases very sharply near both boundaries. When 
the concentration (j) increases, these maxima are quickly 
lowered, and the number density p(r) becomes almost 
uniform. 
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IV. SUMMARY AND DISCUSSION 

In this work we continued our investigation of the en- 
tropic, fluctuation-induced, attraction between adhesion 
bonds in supported biological membranes. The focus 
here was on the possible role played by the membrane 
boundaries as nucleation seeds for the formation of ad- 
hesion domains. Residence of adhesion bonds near the 
cell boundaries is thermodynamically favorable because 
it lowers the entropic cost associated with the suppression 
of the membrane thermal fluctuations around the adhe- 
sion bonds. We started our investigation by looking at 
the pair interaction between an adhesion bond and a cir- 
cular membrane inclusion, and calculating the distribu- 
tion function of the adhesion bond around the inclusion. 
In the case of an open outer cell boundary, the probabil- 
ity density of the adhesion bond decays algebraically with 
the distance from the surface of the inclusion. Although 
the mean pair distance grows linearly with the size of 
the membrane, the distribution function is such that the 
adhesion bond spends most of its time in the immedi- 
ate vicinity of the inclusion. When the outer membrane 
boundary is closed, the probability distribution per unit 
area is symmetric with respect to the mid radius. In 
this case, the adhesion bond is likely to be found near 
the outer boundary, which is much larger than the inner 
one. Diffusion of the adhesion bonds between the two 
boundaries is limited by the existence of a substantial 
free energy barrier in the middle of the membrane. 



Does this analysis imply that the cell boundaries can 
serve as nucleation seeds for adhesion domains? Per- 
haps yes, but probably not due to the fluctuation-induced 
mechanism alone. Our simulation results show that even 
at small densities of adhesion bonds, the entropic attrac- 
tion toward the cell boundaries are entirely self-screened. 
This observation is consistent with our previous study of 
systems with periodic boundary conditions (no bound- 
aries), where it has been found that the fluctuation- 
induced interactions are not sufficient to allow the for- 
mation of adhesion domains. In other words, without 
additional attractive interactions, the problem of aggre- 
gation of adhesion bonds is similar to the problem of gas 
to liquid condensation above the transition temperature. 
Introducing a nucleation seed into such a system without 
changing the temperature is not going to drive the system 
into the condensed phase. We leave for a future work the 
problem whether domain formation occurs at high den- 
sity of inclusions whose inter-distances are smaller than 
the screening length for the adhesion bonds. Another as- 
pect of the problem that needs to be explored is whether 
the existence of direct residual attractive interaction be- 
tween the adhesion bonds and the inclusion can facilitate 
cluster formation, even in the absence of direct attractive 
interactions between the adhesion bonds themselves. 

We thank Haim Diamant and Thomas Weikl for help- 
ful comments. This work was supported by the Israel 
Science foundation (Grant Number 946/08). 
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